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ABSTRACT 


The intonation of a speech is an indication 
of the type of statement of the sentence eg. an 
interrogative, a declarative etc. Many important 
inference can he deduced by knowing the intonation 
pattern in a continuous speech recognition system. 
Intonation are of two types: word intonation and 
sentence intonation. A method has been developed 
to find out the sentence intonation here. 

Vibrations of vocal cords result in the fundar- 

mental frequency (f ) of speech. The variation of f 

^ ' * 

with time is the intonation contour. 

The variation of f has been evaluated with 
respect to time by digital computation. The analog 
speech signal after conversion to digital form is 
processed through a bank of band pass filters simulated 
in the digital computer. The output of this bank is 
stored in a C-matriz. The contents of C-matrix indicate 
the intonation pattern of the utterance. 
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CHAPTER 1 


INTRODUCTION 

1.1 Motivation 

With the vast development in the digital computers, 
the communications with machines have already made the 
transition from science fiction to a practical reality. 
Modern computers now can handle problems considered to 
be non-trivial and complex by human beings. This has 
opened the path for more versatile, ’natural r man-computer 
communication techniques. Thus recent years have seen the 
progress of time sharing systems for rapid man-machine 
interaction* graphical display of pictures, graphs, 
languages on cathode ray tubes; and is discussed herein 
an approach for man-machine communication by speech(voice) . 
Among the various channels available for communication with 
machine, the speech has an edge over the others [23]. 

A few advantages of speech channel are: it 
provides most effortless encoding of all output channels, 
has higher data rate compared to other output channels, 
prefarable for spontaneous output, does not tie up hands, 
eyes, feet or ears; can be used while in motion, requires 
inexpensive and readily available terminal equipment. 

Since the future man-machine communication would 
be desirable through voice, due to the above advantages 
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and since a computer is a most versatile machine we are 
interested here to findout how digital computers can "be 
used in speech recognition systems. 

1.2 Role of Computer in Speech Recognition 

The speech recognition by computers can be viewed 
from two angles. 

(a) The contents of the message 

(b) The linguistic aspects of the message. 

In the former, the words used in the speech are analysed 
to gather certain information for the purposes of under- 
standing them (called the semantics). For example >some 
particular words (say, digits) are used to perform certain 
functions [28 ]. These words spoken by a person could be 
stored suitably in a table or array in a computer memory 

and later on, when spoken again, they could be compared 

the 

with the elements of /table and if matched, certain chain 
of actions can be initiated. The linguistic aspect^ like, 
the stress on vowel, tone and intonation are not utilized 
herein. However, it forms part of type (b). Though much 
work has been done in this aspect of speech recognition 
[1>4#9»21] yet because of the 'uncertain nature* (or 
dependency on number of factors like emotional state, 
sex, age etc.) of the speech, much more is required to 
enhance the man-machine communication. Of immediate 
concern to us here is the intonation. We would venture to 
extract the information about the intonation of speech 
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signal which can then ibe classified to indicate various 
types of statements. 

1.3 Definitions Intonation 

Intonation may be defined as the gross variation of 
the voice fundamental frequency over a sentence or phrase 
where a phrase or sentence is made up of one or more word 
strings in a given language. 

Both the phrase and sentence intonation carry important 
information from speech recognition point of view. What we 
shall be considering in the present work is the sentence 
intonation. As will be seen shortly the various intonation 
patterns would be obtained for different types of statements 
and hence their classification. 

1.4 Speech Production Mechanism 

The generation of intonation of an utterance is 
organized in parts, in terms of certain synchronized patterns 
of muscular activities of the larynx and the respiratory 
system. Intonation is therefore perceived in terms of complete 
contour of fundamental frequency and amplitude. The sounds 
are produced by a few basic mechanisms obtained in various 
ways. A regular periodic vibration is generated through the 
action of the vocal cords, which are two bands of elastic 
tissues in the larynx. They may be opened to permit free 
breathing or brought together to produce various types of 
voice qualities. At a given moment, they are closed tightly 
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with narrow opening between then and if the air behind is 
moderately compressed by the contraction of thoraic cavity, 
the air pressure forces the vocal cords apart. A small 
amount of air passes through them and they again are closed 
because of their elastic nature. And thus the alternate 
closing and opening of the vocal cords give rise to voice 
fundamental frequency(f Q ) . 

Based on the above facts the concept of source-filter 
model of speech production [7] was developed, as shown in 
Figure 1. 

Accoustic Accoustic 

source filters 



e(t) f(t) 


Figure 1.1; A simplified model of speech production 
system. 

The vocal tract is excited by vibrations of the vocal 
cords modulating the air stream passing through the larynx. 
The variation in volume-velocity of the air stream is the 
quasi-per iodic pulse like waveform that is ultimately res- 
ponsible for the fundamental frequency of the voice. The 
vocal tract consists of oral, nasal and pharyngeal ■ cavities. 
It may be considered as a linear filter network with a 
fundamental frequency waveform as its input and the radiated 
speech signal as its output [8], 
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Prom the communication point of view [20] physio- 
logically in general, for males the fundamental frequency 
varies from 70 Hz to 160 Hz; for females 140 Hz to 320 Hz, 
and upto 500 Hz for children. Health factors affecting the 
f here would he assumed to he inconsequential. The effect 
of emotions on the fundamental frequency f Q have been studied 
extensively by others [20-]. 

Sentence intonation is nothing hut the contour , of 
the fundamental frequency f Q . Therefore measurement of f Q 
over the period of utterance would result in the formation 
of intonatuon patterns. 


1.5 Types of Intonation Patterns 

The various types of sentences may he distinguished 
from the study of the variation of f during the utterance. 
Terminal glide in the f Q gives an important clue to the 
type of the sentence. Questions are often distinguished 
from statements hy a terminal rise in f as against a 
terminal fall in a statement. However, a question may 
also he distinguished hy a comparatively high f Q through- 
out the utterance P-6], Spectrographic analysis of speech 
also have shown that, questions tend to he spoken on a 
higher f Q than statements usually ending in a moderate 


rise P-6], These facts suggest that not only the direction 
and range of terminal glide (rise and fall) hut the shape 
and level of the entire contour affect the listner's 
judgement as to the type of sentences perceived. Broadly 
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speaking the sentences could he divided into the following 
major categories s 


1 . Normal statement 

2 . ..Questions 

3 • Anger 

4. Delighted surprise (exclamation). 


As sho-vm elsewhere [12], the average . f Patterns for the 
intonation of above 4 types could be illustrated as in 
Figure 1.2,. The utterance chosen was ’five thousand six 
hundred ten* spoken in different ways. The variation of 
f at the terminal glide is marked in the case of question 
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Figure 1.2* Average fundamental frequency patterns of the 

intonation contours.., . 
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Based on these patterns we can analyse the fre- 
quency component of an utterance with respect to time 
and defending upon the type we may classify the type of 
sentences. 

1 • 6 A Recognition System 

It will he appropriate here to see one example 
of Man-computer communication system through voice and 
orient our problem of intonation in this environment. 

Let us consider a block diagram of a general 
purpose man-computer communication system through voice 

[26] as in Figure 1.3* 

It is clear from the block diagram that the input 
speech after digitization is stored in the computer 
memory. This is then analysed through suitable spectrum 
analyzer. The output of the analyzed speech is recognized 
as a string of sounds (phonemes) and then segmented to form 
sentences and suitable classification is done thereby in 
phase V; the message reconition interpreation of the speech 
is done and appropriate responses are generated suitably 
and outputed as desired bythe interactor. Careful study of 
literature on speech analysis [ 8 ] shows that even with 
different types of spectrum analy z ers often it is diffi- 
cult to resolve the ambiguities of recognition of a diffe- 
rent classes of speech sounds, e.g. a vowel /0/ and a glide 
/w/ might be confused on the basis of their gross spectral 
information - formants, duration etc. To resolve these 
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difficulties, one may think of two of the following 
approaches : 

(1) ask the speaker to repeatedly speak the message 
very slowly and clearly. However, this may not 
always he possible, 

(2) take the cues from other acoustic-linguistic aspects- 
say, intonation contours. 

The second seems to be viable because even though the 
recognition of a particular speech sound may be difficult in 
a given sample, since an intonation type is always present 
in a speech, this can be utilized suitably. To illustrate 
this point, let us take an example: 

Suppose a normal man with normal organs of speech 
(vocal set up) speaks a sentence in English: say (we have 
transcribed in normal orthography here). 

When did you come ? 

If on analysis, we note that the speech segment /w/ 
could not be clearly, and unambiguously (say undiff erentiable 
from /0/ or /v/) identified, what one can do is as follows: 

1. Look at the intonation contour. 

2. Does it raise at the end of the utterance ? if yes, and. if 
the second segment is identified as /h/, and a vowel 
after it, then the utterance begins with wh - (who, 

when etc. ) . 

3. Otherwise call other recognition algorithms. 

Without elaborating further, we note that in English, 
normally a qu'erry sentence ends with raising terminal 
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intonation, irrespective of whether a quefry word, such as 
what, when, which, who, occurs at the beginning or not. 

Thus one can construct a complex algorithm suitably incor- 
porating these facts. The important point to be noted here 
is, that the intonation analysis of a speech input could be 
really used to call appropriate types of spectrum analysis 
techniques (hardware or software wise) and thus enhance the 
recognition procedure enormously. Secondly, while the normal 
speech analysis techniques restrict their attention to the 
recognition of vowels, syllables/words, this approach would 
open procedures for continuous recognition of speech. To 
our knowledge, we are first to introduce this approach in 
speech recognition system. 

In Figure 1 . 4 , we have provided the modified version of 
Figure 1.3 system incorporating the intonation analysis 
procedure before the spectrum analyzer is used. 

In this thesis an attempt has been made to establish 
the feasibility of the analysis, recognition and classifica- 
tion of intonation patterns as a first step in the direction 
of total continuous speech recognition. Using a bank of 
simulated digital filters for this purpose, a preliminary 
investigation has been done with the synthesized signal 
samples. In Chapters 3 and 4 digital filter simulation done 
on IBM 7044 is described and the problems encountered and 
the suggestions for further work are given in Chapter £. 



CHAPTER 2 


REVIEW OP LITERATURE 

Speech research studies have gained momentum 
since I960. Comprehensive survey of literature on the 
automatic speech recognition system exists [ 17] . This 
can he further supplemented suitably to include materials 
published upto 1973 (Appendix A). Most of the studies can 
be broadly classified as follows; 

(a) Isolated spoken vowel recognition 

(b) Syllable recognition 

(c) Humeral (digits) recognition 

(d) Word recognition 

(e) Speaker recognition 

(f) Prosodic feature recognition (stress and tone). 

Apart from the areas of speech recognition 
different techniques, hardware and software have emerged 
in hthe recent past. A few broad systems can be summarized 
as follows; 

-(a) Hardware recognition system 

(b) Hardware cum computer combination 

(c) Adaptive (learning) recognizers 

(d) Usual accoustic analysis 

(e) Simulated recognizers. 
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In the following brief outline of the various 
recognition systems we shall mention only the most 
significant ones and rest of them may be found in the 
literature [ 17]. This however does not mean that we 
underestimate other approaches or overestimate the 
contribution of the systems to be mentioned. 

Before we deal with the various speech recognition 
systems, it will be in order to mention briefly what we 
mean by speech analysis which forms the preliminaries 
required to understand the ensuing survey. l 

The usual way of visual portrayal of speech sounds 
is the use of sonagraph or spectograms. Spectograph (or 
Analyzers) produces a permanent spectogram of a complex 
signal. Developed at the Bell Telephone labs [15] in mid- 
forties, it is used in extracting many speech features such 
as formant structure, voicing, friction, nasality and pitch. 

2 .1 Hardware Systems 

An example of the use of hardware system in the 

recognition of conversational speech is given by Suzuki! 5 0] usii! 

a 26-channel spectrum analyzer covering frequency range from ( 

200 Hz to 5>900 Hz. The channel outputs were separately 

grouped and connected to individual vowel decision circuits. 

The speech was segmented into voiced and unvoiced segments 

and envelope intensity measurements were also used. The 

' ■ ■ . | 
final classification was made by observing the phoneme 

most frequently recognized. 
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A device for recognizing Japanese digits spoken 
in isolation [21] was based on number of voiced interval 
in each utterance, the presence of an unvoiced consonant 
and formant frequencies and frequency slopes at three time 
points in each voiced interval. 

Another promising approach to real time speech 
recognition using analog threshold logic (ATI) gates has 
been developed and the RCA laboratories [22 ] . These gates 
give zero output until sum of excitatory and inhibitory 
inputs is less than a specified threshold. Once this 
condition is satisfied, the output is linearly proportional 
to the sum. The basic ATI gates are used to abstract 
relevant features from the output of the spectrum analyzer. 
Some of the various features extracted are formants, zeros 
positive and negative slopes etc. Such a feature detector, 
instead of a categorical decision as to whether a feature 
is present or absent, gives a quantitative measure of the 
feature. Outputs of these feature detectors are combined 
in appropriate manner to perform recognition. 

A system for recognizing Italian numerals [11 ] uses 
a 17 channel spectrum analyzer followed by threshold 
detectors which reduce the range of outputs from each 
channel to binary statements. These signals are then 
sampled and connected to separate circuits for recognition 
of each digits. 
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2 . 2 Hardware- Computer Combination 

Denes f 6 ] gives a system for recognizing the ten 
digits from complete word patterns. Outputs of a 17 channel 
spectrum analyzer were sampled and recorded on a magnetic 
tape. This is then inputted to a computer to form a time- 
frequency patterns for number of utterances of each word and 
stored as reference. Cross-correlation process was used to 
compare the unknown utterance ’with the stored patterns. Best 
match pattern was chosen as that of the unknown signal. 

Another system due to Gold [ 13 ] uses a 16 channel 
spectrum analyzer, a pitch extractor and a voicing detector. 
The computer program segmented the words into approximate 
phonemic units on the basis of voicing magnitude and spectral 
information. Further spectral durational and intensity 
measurements were then made on a group of five segments 
centered on one segment which was defined as stressed. A 
scoring method, based on the similarity with previous 
measurements on known words was used to classify the word. 

2.3 Adaptive Recognizers 

The system is trained first for an utterance by a 
single talker and then used for recognition purposes. A 
recognition system [ 5 j using adaptive threshold elements 
incorporates a 15 channel spectrum analyzer, a 13 level 
amplitude quantization, and the output is sampled every ten 
m.sec. This is then fed to ten adaptive threshold elements 
simulated by computer. 
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Another system [ 9 ] divides the frequency range of 
300 Hz to 3000 Hz into ten channels, each of which is 
envelope detected. Two more signals are formed from the 
energy above 3 KHz and the overall spectrum. These 12 
signals are sampled for \ sec. after the onset of a word 
and the resulting 192 samples were punched on to a paper 
tape to be fed to computer subsequently. Word recognition 
is by means of linear decision hyperplanes which are adjusted 
to give optimum performance during a training session. 

Recently Clapper [ 2 ] has given a scheme for adaptive 
memory technique to adjust to each speaker. Separately 
spoken individual words can be automatically recognized using 
a two dimension pattern of spectral density vs, a nonlinear 
time base. The pattern for a given word differs from person 
to person and must be adaptively learned by the machine for 
each speaker. Simple circuitry is given that learns a word 
with a single utterance and recognizes it thereafter. 

2.4 Usual Accoustic Analysis 

A highly successful recognition system based on 
direct analysis of speech wave has been developed by 
Reddy [27]* Direct processing of speech wave for segmen- 
tation and feature extraction is the significant aspect 
of this system. Speech wave after normalization is 
subjected to a segmentation procedure which identifies 
sustained and transitional segments. The main parameters 
used for this segmentation are intensity and zero crossing 
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rate of consequetive 10 msec, sections of the signal. 

Minimal segments (10 msec.) are grouped together if they 
are to be accoustically similar, similarity being defined 
on the basis that intensities of adjacent elements do. ; s 
not vary by more than a specified tolerance interval. 

A vowel recognition system [30 ] performs the 
spectral analysis by means of an 'analogue ear*. The 
output of analogue ear was sampled at a high rate in 
4 msec, bursts synchronized with the onset of each glottal 
cycle.- Analysis was performed by a computer by means of 
correlation operations. 

for 

Apart from these studies excepting /one study on 
intonation analysis [10 ] other study of intonation pattern 
exists only from the psychological point of view, for 
example from the perceptual point of view the relation 
between the movement of fundamental frequency f 0 with 
respect to time in giving rise to various intonation 
pattern has been brought out in [12]. This does not involve 
any recognition by a machine but throws much light on the 
perception of the intonation patterns. 

Thus, we realize the necessity or desirability of 
using intonational analysis through a computer system which 
does not exist elsewhere. We describe our approach in this 
regard in the ensuing section ./ 



CHAPTER 5 


STATEMENT OF PROBLEM AND PROPOSED SOLUTION 

It has been observed from the literature survey on' 
the speech recognition techniques that very little signi- 
ficant work has been done in the field of intonation 
analysis. Total continuous speech carries important in- 
formation in the form of intonation contour. Hence the . 
need arises of the analysis, recognition and classification 
of intonation patterns which would be a step towards the 
enhancement of the work on speech recognition. 

The aim of this thesis is to obtain a procedure 
to be implemented on digital computer to analyse and 
recognize the intonation patterns of the specified 
sent ences. 

The variation of the intonation contour is essen- 
tially a plot of fundamental frequency f Q of the speech 
signal with respect to time, f varies during the utterance 
of the sentence, hence its evaluation versus time by digital 
computation constitutes the crux of the problem. It has been 
mentioned in Chapter 1 that f lies in the range of 70 to 
500 Hz. for male, female and children speeches. However, 
for our purposes, we would limit this range from 70 to 300 Hz 
which includes the f Q signals cf male and female speakers 
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In order to process the speech signal in the 
digital computer, it is essential to have the analog 
speech signal converted to digital signal. An A/D 
converter would give the corresponding digital signal 
for speech. The extraction of the information at out f 
at various time intervals can he done h y having a hank 
of hand pass filters, covering the entire range of 70 to 
300 Hz. It would he desirable to have the bandwidth of 
filters as narrow as possible, hence giving rise to a 
large number of filters in the hank. Compromise between 
the bandwidth and the number of filters would have to be 
done depending upon the filters used. We would be using 
the digital filters whose stability is a factor of the 
word length of the computer, the amount of round off error, 
and truncation error and the location of the poles of ■ 
transfer function with respect to the unit circle in 
z domain [25]. Therefore the number of filters and the 
bandwidth chosen would be decided upon by: 

(a) Stability of the filters. 

ft) Time requirement for computation 

(c) Memory requirements of the program. 

The system which would be used for intonation 
analysis-recognition is shown in Figure 3»1*. The 
following assumptions are made in the analysis. 

(a) The .continuous signal to be analysed has the 
explicit boundaries (of the sentence) defined ' already. 
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(To) The duration of the sentence spoken is approximately 
of 2.4 seconds. 
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Figure 3.1: Proposed analysis recognition systen 
for intonation patterns. 

The system comprises of three main blocks: 

(a) Converter 

(b) Analyzer 

(c) Recognizer. 

The converter constitutes a transducer and an 
A/D converter. The transducer produces appropriate 
electrical signal corresponding to input speech, to be 
fed to the A/D converter. The A/D converter gives the 
digital output for the analog (speech) signal input. The 
speech signal would be sampled at chosen regular intervals 
and a binary coded signal would be available corresponding 
to these sampled values. The sampling should be done 
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at a rate twice the maximum frequency present in the 
speech. The output of the. converter is stored on a 
magnetic tape to "be processed on the computer later on. 

The blocks analyser and recognizer are simulated 
within the computer. The digitized speech is first passed 
through the low pass filter. 'A' to give signal sequence' 
containing frequencies below 500 Hz. This sequence is 
then passed through the band pass filters BPF 1 to BPP^ 
successively. The output of the filter bank is stored 
in the C matrix (code matrix), which represents the time- 
frequency relation of the speech signal. This is the 
intonation pattern of the utterance or the speech signal. 

A new set of computational procedure is required by 
which the contents of C-matrix can be classified by the 
recognizer. The standard intonation patterns for various 
types of sentences are stored in the computer memory. By 
selecting suitable algorithms the C matrix contents should 
be compared with the stored pattern. Decision on the type 
of the pattern should be made based on the above comparison. 
This aspect of the system requires a further investigation. 

However, in this thesis we are restricting ourselves 
only to the analyzer part and rest of the system is assumed. 
In the next chapter, we give the details of a simulated 
analyzer. 



CHAPTER 4 


SIMULATED ANALYZER 

Processing of the digitized speech can best be 
analysed IT or various frequencies by passing the signal 
through a bank of digital filters simulated in the 
computer. What follows in this chapter is a discussion 
on realization forms of digital filter transfer function, 
the design criterion, simulation of filter bank and the 
choice of these filters for our purposes. 

The term digital filter refers to a system which 
executes an algorithm by which a sampled signal or 
sequence of numbers x(nT), acting as an input is trans- 
formed into a second sequence of numbers y(nT) termed 
the output signal. 

One such algorithm used very often is 
N M 

y(nT) = a.x(nT-iT) - *> ' b^yCnT-iT) (l) 

1=0 1 i=i 

Knowing the present and (N-l) past values of input 
sequence x(nT) and the M past values of output sequence 
y(nl) the present value of the output can be computed. The 
difference equation (1) in time domain gives the transfer 
function of the digital filter in z domain as 
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JL .1 

£ a i * 1 

H(z ) = — - (2) 

M 

1+ b.z"* 1 

i=i 1 

where z” 1 represents the unit delay operator. 

The transfer function H(z) given hy (2) leads to 
two types of filter realizations: (a) where at least one 
a ± ■ and one b^ is nonzero - gives recursive digital filter, 
and (b) where all tu are zero, which gives nonrecursive or 
transversal digital filter. 

Recursive digital filter output depends not only 
on the input sequence but also on the previous values of 
the output. Whereas the transversal filter output depends 
entirely on the input sequence. Both these types have 
distinct characteristics. The transversal filter possesses 
excellent phase characteristics, but it requires very large 
number of terms (hence more number of multiplication and add 
operation to obtain a sharper cut off. However the recursive 
filter requires few terms to give a much sharper cut off. 
Thus to obtain same cut off recursive filter requires less 
computational efforts. 

The required transfer function H(z) for digital 
filter can be obtained in one of the wtwo ways: (a) obtain 
the transfer function in s-domain for the corresponding 
analog filter and apply suitable transf ormation to obtain 
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z-domain transfer function, []b3p4+ (b) obtain the 

transfer function of the digital filter directly in 
z-domain [25 ]• Both methods have their advantages. In the 
first case, the use could be made of the extensive work 
done in analog filter synthesis. Among the various choice 
of filters, a few are Butterworth, Chebychev and elliptic 
filters. All these filter characteristics are approxima- 
tions to a desired rectangular pass band. The Butterworth 
kind achieves this through monotonic amplitude versus 
frequency characteristics. By allowing a ripple in pass- 
band the Chechychev kind using the same number of poles 
and zeros can achieve sharper:-. cut off* Elliptic filters yield 
even Sharper cutoff than Chebychev for same network complexity. 


Synthesis techn: ; :;ies for analog filters yield transfer 


function in s-d amain: 


K 




H(s) 


k=o 


H dnS- 
11*0 1 


E. < 1 


( 3 ) 


The transfer functionoof digital filter can be 
obtained from this by applying one of the following 
transformations : 

(a) Standard z-transf ormation 

(b) Bilinear z-transf ormation. 

In first case transfer function H(z) is evaluated 

log z 

using the transformation s if— » where T is sampling 

frequency. 
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T&us 


H(z) = H(s) 


>s= 


log e z 


(4) 


This is conveniently done by breaking H(s) into 
partial fraction and then making the substitution 


1 

s+a 


T 


n -1 -aT 
1-z e 


Standard z transformation is good in preserving 
the impulse response but its frequency response suffers 
from the aliasing error. This error is due to the fact 
that H(s) is not band limited. However, if.it is band 
limited, i.e. H(s) = 0 for, wj > w g /2 ; w g being angular 
sampling frequency, then aliasing error would be absent. 

This aliasing error, can however be eliminated by 
using the bilinear z transformation or z-form. This 
requires the substitution 


1-z 


-1 


1+z 

This transformation, eliminating the aliasing 

error, carries a nonlinear warping of the frequency scale 

due to the fact that the whole of the frequency scale axis 

* 

of s plane is compressed into a portion j w j <T w s /2. This 
can be compensated by prewarping the critical frequencies 
of the analog filter in such a way that application of 
z-forms win, sfeift these frequencies back to their original 
values • 
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Thus 

D. 

^ = tan ~~ (5) 

gives the analog frequency for corresponding digital 
frequencies wu [25]. 

-*■' jL 

Another transformation matched z transformation 
could also be used for obtaining H(z) [14]. 

It is ufeual practice to design a low pass filter 
first and then using suitable transf ormat ions , high pass, 
band pass or band stop filters could be obtained. These 
transformations could be applied either in analog filter 
[3l] or in digital filter [ 3 ] transfer function. 

It has been shown [ 3 ] that if H(z) is the transfer 
function of low pass filter, then H(l/f(z)), where 

-e f \ 2 ( z —a ) ■ 

■ (1-a z ) 

acquires a band pass characteristics. 

Here a is so chosen as to give appropriate cut off 
of the band pass filter, a is given by 

cos k(w ? + wo ) 

a - cos(2k w ) = (6) { 

u cos k(w 2 - w-^) 

where w Q = centre frequency of band pass filter 
w 2 = upper cut off frequencyccy.BP filter 
w^ = lower cut off frequency hflBP filter 
k = pr ji\ = angular sampling frequency. 
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The cut off frequency for the corresponding low pass 
filter is given by 

w 0 = W 2 - W 1 

Realization of Filters 

A recursive filter represented by (l) can be 
realized (or simulated on computer) using either direct, 
cascade or parallel form as shown in Figure 4*1. 

The direct form is the direct realization of the 
expression of transfer function nwhereas cascaded form 
is obtained by expressing transfer function as 

H(z) = (z ■*•) x (z ■*") x ... H^(z 1 ) 

where each of the subfilters includes a subset of the 
poles and zeros of H(z~"^). Parallel form is obtained by 
expressing H(z”" L ) in its partial fractions. 

Direct form requires a larger amount of memory and 
greater accuracy in filter parameter determination. This 
is seldom preferred to the cascade or parallel realizations. 
However, cascade form is preferred over parallel form 
because the errors due to finite word length in computer 
is pronounced in parallel form. 

System Design 

The speech output of the A/D converter has been put 
on the magnetic tape, for processing later on the computer. 



Direct form H(z ) = — 

i+fCz”- 1 -) 



Parallel form Riz* 1 ) = H^z” 1 ) + H^z"' 1 ) + ... + 


Figure 4.1; Various realizations of digital filter 
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The sampling rate of the signal is 20 KHz using a word 
length of 12 bits. Offset binary was used in coding. 

As has been pointed out earlier, the intonation 
contour is the variation of fundamental frequency ..:f , 
of speech with time. It has also been shown that this 
fundamental frequency f of speech lies in the range of 
70 to 500 Hz. 70 Hz to 160 Hz for male speech, upto 
300 Hz for female speech and upto 500 Hz for children. 
Therefore, in order to evaluate the intonation contour, 
it becomes necessary fhattfehec signal be passed through a 
low pass filter with cut off at 500 Hz. before processing. 
The resulting signal could then be analysed for intonation 
contours. Specifications for this filter would be 
Cut off frequency 500 Hz 

Pass band tolerance 25 percent 

Stop band tolerance 15 percent 

Transition ratio 0.8 

Sampling frequency 20 KHz. 

The type of filter chosen for this is elliptic 
filter. Since it would give the required filtering with 
minimum number of cascaded sections. A computer subroutine 
PROXP due to Sinha [ 30] was used to carry out this design. 
This subroutine returns the coefficients of the digital 
filter which is simulated as subprogram DIG-PIL to give the 
required output. 
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Filter Bank 

The fundamental frequency f of 70 to 300 Hz has 
been divided into contiguous channels of equal bandwidth 
of 25 Hz each giving the 10 channels of following bands. 


70 

- 95- : 

Hz 

195 - 

220 

Hz 

95 

- 120 

Hz 

220 - 

245 

Hz 

120 

- 145 

Hz 

245 - 

270 

Hz 

145 

- 170 

Hz 

270 - 

295 

Hz 

170 

- 195 

Hz 

295 - 

320 

Hz 


Here again the type of low pass filter used was 
elliptic filter and the subroutine PROLP was used to obtain 
low pass filter coefficients. A prototype low pass filter 
was designed first with following parameters: 


Cut off frequency 
Pass band tolerance 
3top band tolerance 
Transition ratio 
Sampling frequency 


25 Hz 

20 percent 
10 percent 
0.9 
1 KHz 


Sampling frequency in this has been chosen as 1 KHz as the 
maximum frequency present in the signal is 300 Hz. Achieving 
this, however, is no problem since 1 in 20 samples of the 
input signal can be taken to suit this condition. Transfer 
function of the prototype low pass filter is obtained as 


U” 1 )- 


n 


a il +a i2 Z + a i3 Z 
1 b il +b i2 z 1+b i3 Z 


where n is the number of cascaded section. 


( 8 ) 
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For further discussions let us take n = 1 , which 
however can he extended to cases n 1 without difficulty. 

Convex-ting this low pass filter to hand pass filter 

1 - z ~ 1 ( '■ 1 \ 

hy the transformation z~ J " — * — ZSLI j_ n (3) 

(1 - az -1 ) 


a u^i 2 (- )+ a 13 t- 

, -l x 112 (1 - «z d ) 13 (1 - az X ) , Q , 

H(z - 1 ) = e-. rr~ 3 l =T— TT 0 (9) 

Vi 2 (- *„(- l- lS- = p ) 2 

11 12 (1- a z 1 ) 13 (1 - az -1 ) 


Upon simplification, one gets 

-1 -2 -3 -4 

-t C-, +C-, 0 Z +C-, *Z + C-,„Z " / +C 1 r- z ■ 

tt^_-1n 11 12 13 14 15 

hlz ; = 

d ll +d 12 z ' 1+d 13 z " 2+d 14 z ‘ 3+a 15 a ' 4 


( 10 ) 


where the coefficients in (10) are related to that of (8) 
hy following relations 


11 = 

a ll 

d ll 

12 = 

a(a 12 -2a i:L ) 

d 12 

13 

a ^ a ll~ a 12 +a 13^~ a 12 

d 13 

14 = 

a(a 12 -2a 13 ) 

d 14 

15 = 

a 13 

d 15 


a(bi 2 -2b 11 ) 


u 12 


<x(bi 2 -2b 13 ) 


( 11 ) 


Thus one comes to conclusion after observing ( 11 ) 
that the coefficients of various hand pass filters can he 
obtained from a low pass prototype filter coefficients, 
the only parameter changing is a, which is decided upon 
hy the cut off frequencies of the corresponding filters. 
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transfer function in (10) can be represented 
pictorially as shown below. 



This has been simulated on digital computer using 
the subprogram DIJfIL. The unit delay z ^ represents 
storage for one data sample. 

Testing of Filters 

Before using the filters for actual processing, 
they are tested for (a) Impulse response (b) frequency 
response. 

Impulse response indicates the stability of the 
filter. This could be obtained by getting the output of the 
filter for an input sequence of 1,6,0 ... 0. for a stable 
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filter the output decays in transient state and becomes 
zero in steady state. 


Frequency response of the filter is obtained by 
i * ~*1 , — j w r 

replacing z by e in the expression for transfer 

function, i.e. from (10) 


H(e“ J ‘ wT ) 


c +c (e“^ wT 
11 c 12 ie 


< ^11 +C ^12 ^ 


Up -J2wT -j3wT -d4wT 

_/ +c l*> e +c i 4 e +c 15 e 

3wT^ ^ 3wT 0 -J 4 vT 


+d 13 e 


+c 14 6 


’ +d 15 6 


( 12 ) 


dividing numerator and denominator by e ""^ 2wT and simpli- 
fying gives 
H(e~j wT ) 

[ ( c 2j_ +c p5 )cos2wT+(c 12 +c 1 ^)coswT+c 1 -K 
E (dfi+di^ ) cos2wT+ ( di2 + di4 ) coswT+d-^]{ 

+j[ (c 12 -c lt5 )sin2uT+(c 12 -c 14 )sinwT ] 
+j[ (d 12 -d 15 )sin2wT+(d 12 -d 14 _)sinwr ] 


H(w) for various values of w can be plotted to find whether ! 
the response is the desired one. 

Program Implementation j 

' 

The aperiodic nature of the speech wave necessitates jj 
the creation of a data window. This subset of data (i.e. a 
portion of speech signal) is passed through the bank filters 

respectively number of times till the filter output stabilises . 

' . . ■ [ 

It has been found by experimentation that 5 such repetitions | 

give the desired results. The data window chosen was of 
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50 samples. Once the corresponding output is obtained 
the r.m.s. value of these 50 samples is taken and compared 
with the r.m.s. value of the corresponding input samples. 

To ascertain whether the output of a particular filter 
is present or not, use is made of the fact that filter would 

attenuate the signals outside its pass band. Whereas the 
signals in the pass band would give larger outputs* Thus if 
the ratio of the r.m.s. values of the input and output of the 
filter is above certain threshold then the presence of the 
particular frequency is indicated in that particular subset 
of speech signal. 

Thus the decision taken at the end of the processing 
of 50 samples, based on the value ratio will be indicated as 
1 or 0. If the signal is above the threshold, indicating the 
presence of frequency in that period, would be indicated by 
1, and the absence by 0. 

The next 50 samples are taken subsequently and 
processed similarly, till all the samples are processed. 

The signal is then passed through the subsequent filters* 

The decisions 1 or 0 at various time instants are 
stored in c-matrix. 

The computer program organization and the flow 
charts for the above program are given in Appendix B. 



CHAPTER 5 


RESULTS OP EXPERIMENTS 


The low pass filter A (figure 5 .1 ) for speech 
signal and the hank of filters were designed and simulated 
on the computer. The coefficients (i.e. a^'s and tu’ s ) 
of the low pass filter A and that of prototype low pass 
filter are shown in Table 5.1* The coefficients of the 
various hand pass filters were obtained from the low pass 
prototype filter using equation (11 ) of Chapter 4* These 
coefficients (c^'s and d^'s) are tabulated in Table 5*2. 

The testing of the simulated analyzer was carried 
out in the following manner. 

A signal comprising of five sinusoids of different . 

frequencies was fed to the program package. Pive filters 
in whose pass bands these sinusoids are covered indicate 
their presence by 1 in the relevant row. Whereas the 
absence is indicated by a 0 in the row for which no _ 
sinusoid existed. The results are shown in R.i. Another 
similar signal, but the frequencies of sinusoids staggered 
from that of the first signal was also fed. This results for 
this are shown in R.2. These results establish the efficacy 
of the filter bank. 

Because of the unavailability of the digitized 
signals were synthesised by mixing several 


9 
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sinusoids together so as to resemble the spectrum of 
the speech signal. Sinusoids of different frequencies 
appeared at various tine intervals in such a signal. 

Two such signals were synthesised and were analyzed by 
the program successively. Results for these are shown 
in 1.3 and R.4 displaying the c matrix. They indicate 
the presence of frequencies in the various time intervals 
and in corresponding bands as per the synthesised signal. 

It therefore can be inferred that similar type 
of results would be obtained when a digitized speech 
signal is given as input to this program package. 









































































fable 5.2 (continued) 
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CHAPTER 6 


CONCLUSION 


The intonation contour is an important feature 
of speech. Its evaluation would prove to be a signifi- 
cant tool in the continuous speech recognition systems. 
In the present work a procedure has been developed to 
find out the intonation contour of an utterance and the 
feasibility has been established. Some of the specific 
contributions of the thesis are: 

(a) Simulation of bank of digital filters. 

(b) Use of frequency transformation techniques 

to obtain appropriate band pass filters 
from low pass digital filter transfer 
functions, for the purpose of speech 
spectrum analysis. 

(c) A procedure to evaluate the intonation 
contour of a continuous speech signal. 

The program package developed, gave satisfactory 
results with the synthesised speech signals. However, 
it could not be tested for actual speech signal. The 
magnetic tape containing the output of A/D converter 
(by courtsey of T. I.P.R.Bombay ) could not be used 
probably because of the incompatibility of the systems. 
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Attempts were made to design "band pass filters 
with a narrow bond width (10 Ha). This, however, resulted 
in oscillations at the filter output. It has been found 
that at narrow bands, the poles of the transfer function 
become very close to the unit circle. Due to rounding off 
error in computation, it is possible that the poles fall 
within the unit circle, leading to instability. This 
difficulty was, however, overcome by choosing a BW of 
25 Hz . 

Scope for further work 

In order to obtain better intonation patterns, 
narrow band filters should be used. It is therefore 
desirable to obtain stable narrow band filters. This 
particular aspect requires further investigation. 

Another aspect which requires further investi- 
gation is the selection of suitable algorithms for 
matching the stored patterns with that of the analysed 
ones. Once this has been done, the classification of 
the utterance on the basis of the intonation patterns 
could be achieved. 


Speech Recognition and Analysis Review 
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APPENDIX B 


MAIN 

PROGRAM 


SUBROUTINE 

1CPAS 


FUNCTION 

.pjjzm 


SUBROUTINE 
PRO IP 


.[FUNCTION 1 
I XMAX 


[SUBROUTINE 

; POZE 


FUNCTION 

DIGFIL 


SUBROUTINE! 
QUAD J 


COMPUTER PROGRAM ORGANIZATION 


FUNCTION DIGFIL - 
FUNCTION DIJFIL - 
SUBROUTINE LOPAS- 

SUBROUTINE PROLP - 
SUBROUTINE PCZE - 


SUBROUTINE QUAD - 

FUNCTION XMAX - 


computes the low pass filtered sequence, 
computes the hand pass filtered sequence 
performs low pass filtering (for 500 Hz)on 
the speech input 

computes the coefficients of low pass filter 
computes the poles and zeros of the 
transfer function 

computes the quadratic expressions for 
transfer function 

computes the maximum value of an array. 









PLOW CHART cF COMPUTER PROGRAM 














THIS PACKAGE FIMOS THE INTONATION PATTERN OF AN UTTERANCE 3F D UR/ I 
C 2 • k SEC* 

u VECTOR SLIII) CONTAINS SPEECH INPUT SAMPLES. 

DIMENSION ACS, 2, 3) ,c (6) ,SlI (2 2»0 ) , SSL I < 5,; ) , SLO ( 5 .O , VALR C 4 ) 

, ,VIN (S,4) ,CNT(6) »FR( 5 ) ,8(6 ,2, 5) 

INTEGER DfcCN ( i 0, 40 ) 

COMMON VI N 
NT A 3 E = V 
PI=3. 14159 

C READ SPEECH SAMPLES I NT C SLIII ) r ROM TAP,-:. 

READ { NT AP:.:. ) SLI 

C CALL L3PAS TO ELIMINATE FREQUENCIES ABOVE 5C_ . HZ. 

CALL LOPASt SLI ) 

C DESIGN PROTOTYPE LOWPASS FILT.. R 
i, T R = T R A N S I T I 0 M R A T I U . 

L PB1=PaSS B A i‘i D Tj L K A ij C c . 

L S8T -STOP SAND TOLc RANGE. 

C FSAK=.SAMPLI NS FREQUENCY. 

C' :UTFFL=CUT OFF FREQUENCY OF PROTOTYPE LOWPAS FILTtRC THIS £Q, 

C BAND WIDTH OF THE 2P FILTER) 

DATA TR,PBT,SBT,FSAM,CUTFFL/.9*.2».1,;.000.*25./ : 

P I S M = P I / f SA F 

C GLNcRATS FRrl.yUE.NCl ES FOR BP FI LI R ; 

FRO. 3=7. . 

DO 68 1 = 2, 2 ; 

FR(I )=FR( I-D + CUTFFL 

«S CONTINUE i 

... SUBROUTINE PROLP RETURNS COEFFCI 1NTS OF IP FILTER IN ARRAY AC I, J,k| 
C I -INDEX OF SECTION CASCADE { 

U J-IS .. FOR NUMERATOR i 

C J-IS 2 FOR DENOMINATOR j 

C K-t.2»3 DESIGNATE COEFFICIENTS OF DELAYS j 

L NSEC-NUMBLR OF SECTIONS IN CASCADE . j 

CALL PROLP (CUTFF L, TR, PBt, S8T»PSAM, E» A, MSEC ) 1 

K = i j 

U COMPUTATION BY FILTER-BANK 

C ND. OF F I L T c R S = ?.. 0 : 

DO? CNF I L-i , 1 

C FL- LOWER CUTOFF OF BP FILT.R 

C rH-JPPER CUTOFF OF BP FILTER | 

FH=FR{ NFIL+1) f 

FL = FR( NFI L ) 

ALP= ( COS ( PI SM* ( FH+FL ) ))/COS (PiSM*(FH-FL ) ) 
t COMPUTE COEFFICIENTS BU*J,K) OF BP FILTER, AFTER APPLYING FREQUENCE 

t TRANSFORMATIONS TO LP FILTER. ] 

D025I=i,NSEC | 

D02 5 J = i, , 2 
B(l, J,U=A(I,J,i» 

B( 1, J,2)=ALP*{A(I, J,2)-2.*A (1,0,1)) 


t 



B( I, J, ) = {A (I, J, A )-A{ I, J*2) +Ai I, J * 3 ) )*ALP**; -A( I, J» 
B ( 1 » J » 4 ) =ALP*( A( I j J? l ^ A (1 5 J ? S ) ) 

B(It Jt5}-A(I,J,3) 

....5 CONTINUE 

DOB 5 1 = 1 j NSEC 
D0i>3 J = i » 2 
D035 <=2,5 

B ( I , J » K ) =B ( I , J i K ) / o ( . *; J s ) 

CiNT ( I ) =B { I ,1,* ) / B i i i ,, j ) 

.,5 CQNT I M U L 

D036 I=i ,N$EC 
DOSS J= i , 

8 ( I j J , » ) = I • 

X. 6 CONTI MU: 

C INITIAL IS,. VIN(1,J)= . 

D371 I=; ,6 
D37iJ=; , ' 

VI M( I , J ) = .0 
Vi CONTINUE 

MB = l 
N 0 =. 

DOTE I W M 0 W = 1 , 2- l 

C READ 5 j SAMPLES OF INPUT INTO SG, I U ) . 

00730= ,5 ■ 

SGLI ( J ) = SL I ( MB ) 

NS=N B+ 

V3 CONTINUE 

Dim IRLP=, ,S 

c OUTPUT OF BP FILTER IS IN ARRAY SLO(I). 

DcmisAM=:. , 5 
X = SG LI ( ISAM) 

SL0( ISaM)=DIJFIL{X,NS£C» 6 »E ) 

'.4 CONTINUE 

C COMPUTE RMS VALUE OF INPUT AND OUTPUT SEQUENCES. 

SUM I N = . 

SUMOT = 0.0 
D095 I SQ=i » 50 

SUMI N= SUMI N+S3LI ( IS0)**2 
SUMOT=SUMOT+SLO(ISO) **2 
95 CONTINUE 

C cVALUTE RATIO OF THE RMS VALUES OF THE TWO SEQUENCES. 
RATI 0= SORT (SUMOT/ SUMI N) 

VALR ( I WNDW ) = R A T 1 0 

L COMPARE RATIO WITH THRESHOLD VALUE THRS. 
IF(RATIO.GT.THRS) GO TO 9 4 
DECN ( MFIL , I WNDW+1 ) = ’» 

GOT072 

94 DECN(NFIL»IWNDW+l)-i 

» L CONTINUE 

PRINT 134, ( VALR ( IK) , XK=X,38 ) 



;0 CONTINUE 

OUTPUT THE C-MATR1K* 

■ J-T1MU COUNTER? ONE INCREMENT REPRESENTS 

' I- BP FILTER COUNTER. 

PRINT 217 
PRINTS,. .9 
PRINT?.::. 6 
DQ92IM=l ,iO 
K=li -I M 
I L=K 

IH-K+L .... , _ 

PRINTi j'i»FR(IL ) ? FR { IH) > ( DEC N{ J ) » J = i J 

92 CONTINUE 
PRINT' J 

prints.:: 

,02 FORMAT! XjlGti, 3.3 i 

FORMAT! X»-.9F3. 2) 

018 FORMAT l X f ^FILTER BO 10*/ ) 

FORMAT! / »63X,*TIME *») 

03 FORMAT { L HO ? 2 P 3 . 0 » 2 X » 3 9 1 3 ) 

,.. .9 FORMAT ( , (/)) 

/2- FORMAT t ///» 25X » * R *v C.-M A T R I X FOR 

2PATTi-RN*) 

:- l '5 FORMAT { X » _9F6.Z?4Xt«iP&.i I 
227 FORMAT O. HI. ) 

STOP 

END 


53 MSEC. FOR KHZ. 


k N I N T 3 N A T I 



*I3FTC 

FUN: 1 1 DM D 1 JFI L { X } N, fi , E ) 

C COMPUTES 'i He BP FILTERED OUTPUT OF THE INPUT SEQUENCE. 

0 VINCI » J ) - TEMPORARY STD RASE FOR cQUENCE 
L J- REPRESENTS DELAY NUMBER 

DIMENSION 8 ( 6» 2 » 5 ) ,t(6) ,VIN(6,4) 

COMMON VIM 
0 Q i I = L f N 

P= C X-3 C I , 2 , 2) * VI N C I , . ) - C ( 1 , 2, ■- )*V I N< I , 2 ) -B C I * 2 , 4) *VI N l I , 
iBCl»2,5)*VIN(I,'+> )/E!I,0. ) 

X= ( 3 ( I 1 i . , ) *p+ B U , t ) * V I N ( I , ) + B { ! , l , B ) *V I N C 1 , 2 J + 6 (I , , 

2 ) .+ B ( I » : f 3 ) * V I :H J»4))*EU )* 2 ... n 6 92 / D . $ 27 A 9 7 

vinci, n)=vin(I f ; ) 

VI N( 1 ,'■)= V I N ( I ,2) 

VIN{ I »*. :)= V I N ( 1 } ) 

Vinci,; )=P 
* CONTINUE 
U I JF1L=X 
RETJRN 
END 


£ ) - 

s~)*VIN 


I 



SUBR 3JT I !ME PR3LP (DM ,TR,R;»R ,OMS , E« A » MI ) 

COMPUTES THE COEFFICIENTS OF THE CASCADED TRANSFER FUNCTION. 
DIMENSION L( 6 ) » A { 6 , 2 , 3 ) » G AHA { 21 ) , ZERO ( 2'. ) , UP ( E, ) , K P { i . ) 
COM 5 LE K POLE! 2 ) ,GUM 
EP=S3RT( IR1*Z* )-<R .-MU 

AL = S QR T { „ . -RL*R Z ) / i E P*R c ) 

0M2= O'i /TR 
XK = l ./(?;. #OMS ) 

C3N= T AN { 6.28- L L *XKKM. ) 

DON= a . /CON 
GON= D3N*DO:i 
3A = 4 L 
1 = 0 
1=1* i- 

GAMA (I ) = GA 

G A ~ ( + , * ) / ( A • ^ Q R ‘ ( ) 1 

R 0 M i = ( 0 n $ / 3 . i 5 9) * A T A N { C 0 N* G A ) 

I F ( R 0 M . GT . 0 M c ) 3 DT 0 ■ 

I R = 0 M 1 /ROME 
NT={ 1**1 )/2 

CALL PGZE (EP, GAMA, i, ZERO, UP* P) 

CALu QUAD (DON, CON, NT, ZERO, UP, \P,e,A> 

J = 

Xc=XMAX( l ,NT) 

J = J +■ ‘i 

IF( = ( J ) .-Q.XUSOTQ ; 

SOTO 2 

£ ( J )==( J]*(i.-R;. ) 

RSTJRM 

END 



*iBFTC 

SUB * OUT I POZ t ( E P » G AM A , NT * ZE P.Oi U P * XP i 

C COMPUTES THE POLES AND ZEROS OF T He TRANSFER FUNCTION. 
COMPLEX POLE (20) * GUM 

DIMENSION GAMA ( £ i 3 » Z : R3 t £ 0 ) , UP I C G ) » X P ( 2 ) 

X= I. /EP 

ZERO ( i 3 =SQRT ( ( ,+GAM il) i/Z.J 
POLE { t ) -CMP LX( D* » X 3 
X= ZE R3 t i, 3 
GUM=CMPLX ( X» .) 

X=GAMA(i 3 

POLE £ 1 3 =GUM#CS CRT ( £ C M P L . { . ,3 . )+POLtC . 3 ) / ( CMPLXt X» 0. } +P3L 
1 F t NT. cQ. 1 ) GOT J b 
DO 1 1=1, HT 
N J i 3 = £ £**1 ) / : 

X. = G4M4 { I ) 

GUM=CMPLX (X j .3 
D01 J = , , ,-JiC- 
NG=M JI G+J 

ZE R3 { N 3 ) = SQRT { ( { . + £3/2. )*t I. ■- ZERO £ J 3 ) /( X-Z:.RO( J ) 3 ) 

ZER3 ( J 3 = SORT ( £ ( .+ 3 /2. i*C .+ ZLRO { J ) ) /(X + Z£RO( J 3 ) 1 

P3LE ( NS ) =CSGRT ( (CMPL U . , . ) + GUM 3 *t CMPLX ( * . , • )-PQL£( J ) ) 

, /(CMPLXt . ,0.)*(GUM-POLE*J) 333 
POLE ( J ) = C S Q R 7 ( ( CMPLXt L. ,0. ) +GUM ) * t CMP LX { i . , D . ) + POLfe ( J 3 ) 
i/l CMPL X t £.» C .)*{ GUM+ POLE ( J J ) ) 3 
CONTINUE 
N J 1 3 = ( tL & N T 3/2 
DO 2 J = , ni JIG 
UP ( J ) = Re A L t POL L t J ) ) 

XP ( J 3 = A I MAG (POLL ( J 3 3 

CONTINUE 

RETURN 

END 



* IB FTC 

SU8R0J T I Ht QUAD { DON, GOM , NT , ZERO, UP, XP , E , A } 
v COMPUTES THE QUADRATIC tXPRESSIOi-J FOR FILTER TRANSFER FUNCTION 
DIMENSION ZERO (Z 0) , UPt 20 ) *X P( 0 ) , c ( 6 ) , A < 6, 2, 3 ) 

DO I=„,NT 
A{ I,L,i)=GON+ZERO{ I ) 

A(I, 1,2)*- 2 .*< GON-Zi. RO d ) ** 2) 

Ad, i»E) =AtI,L , ) 

YY=JP( 1 ) **2+XP ( I ) 

ZZ=2.*ABStXP(I ) )*D- . 

A{ 1,2,3) = YY-ZZ + GON 
At I, 2, 1 ) = L . * { Y Y - GO : ) 

At I, A, ; ) = YY +ZZ + SON 
S d ) = ( A t I , 2 , », } +A t I , 2 , E ) + A d ,2 , 3 ) ) 

/ ( A ( I , , ) + A { I s I i , , ». ) ) 

CONTINUE 

RcTJRM 

END 



FUNCTION XMAX ( CQL, N ) 

DIMENSION COL { 6 ) 

XMAX=D. 

DOH 1=1, N 

IF((A3S(XMAX)-ABS(C3UI ) J J.GE. . } GO TO 
XMAX=ABS ( COL ( I ) ) 

CONTINUE 

RETURN 

ENID 



*1 B FTC 

SUBROUTINE LOP 4$ ( S LI ) 

o subroutine: computes the lowpass filtered output 

DIMENSION SU(/ L O,C{5»2f3 ), (6) »RIN(6*2) 
COMMON RIJN 

DATA TRL, PBTL,S3TL,FSAML»CUTFFL/.8,.25,.i&,2 
DATA RIN/.2* ./ 

CALL P ROLF ( CUTFFL » TRL * P BTL » SBTL,FSAML,E*C, N) 

DO i l*i ,£ y-)V 

X=SL I ( I ) 

SLI { I) =DIGFIL( \,N, i C 3 

CONTINUE 

RETURN 

END 




kJ kj J* ) 5 u %) 0 / 



*18 FTC 

FUNCTION QIGFI L ( V * .. 1 4 ) 

DIMENSION A (5,2,3) »t(6) ,RIN(6,2J 
COMMON R1 N 
DU I=J*N 

P=(i//MIk,1))-[A{I,,,E)/A(1, , . ) )*RIN( I, 
i - ( A( 1 , 2 1 3 ) / A ( I , I , i ) ) *RI Nt I * 2 ) 

V=E( I )*{A(I» j 1 )*P+A( It ' ,Z)*RIN(I,1 ) 
i + A { I , L , 5 ) * R I N ( 1,2) )*:..! 9X692/ .52 7197.,. 0 
RI NT I» 2 ) =RI N( I » .... ) 

RIN( It t ) = P 

CONTINUE 

DISF1l=V 

RETJRN 

END 



C THIS PROGRAM READS THE TAPE 
C TAPE T3 BE MOUNTED ON 8-4 
C ARRAY ADOUT CONTAINS ONE RECORD LENGTH 
C ARRAY DAT STORES EACH OF 12 BITS. 

C ARRAY PAT CONTAINS DECIMAL EQUIVALENT OF THE! 2 BITS CORRESPOND IN 
C ONE SAMPLE 

€ IRCORD=NO. OF RECORDS TO 8.-. READ 

DIMENSION PAT(4Iv ) 

INTEGER ADOUT ( 1333) ,DAT(4 u) ?A»B,C»D,E 

I RCDRD = 1 

D020IM=1 , IRCORD 

D01I=1 ,1333 

ADOUTI I ) =C .0 

E=2**12 

C READ THE TAPE 

CALL READ ( ADOUT ) 

I FLAG= 

J= y 

D010 1=1,1333 

C SEPERATE THE 36 BITS INTO _2 BITS EACH ! 

A = ADOUT ( I 3 

IF(A.GE.O) G0T05 ' 

I FLAG*! ! 

A=-A 

5 C = A/E : 

8 = A-C*E ; 

D =C / £ 

C=C-D*r 1 

IF ( I FLAG .NE. . ) D=D+E/2 

DAT ( J ) =D ! 

DAT( J+i)=C i 

DAT( J+2 3 = B | 

J=J+3 i 

IFLAG=0 I 

. 0 CONTINUE | 

REWIND 1 ' | 

READ (13 DAT I 

C CONVERSION TO DECIMAL f 

DO 1 5K *1,400 D 
PAT ( K) =DAT ( K) 

PAT( K 3 =10 »* ( PAT ( K 3 /4 96 . ) 

,5 CONTINUE 

PRINT 103, PAT 
23 CONTINUE 

FORMAT (SOX, 17/3 

101 FORMAT { 5X , 1CD12 ) ' 

a 33 FORMAT ( 5X, 20F6. 2 ) 

STOP 

END 



*IBMAP COM V 

C MAP SUBROUTINE TO READ THE TAPE 



ENTRY 

READ 

READ 

SAVE 

A f 2 1 4 


CLA 

3,4 


STA 

10 


ENB 

=0 


TCOB 

* 


RDS 

TB4 


RCHB 

10 


TRCB 

ERR 


EN3 

*0000.- 


RETURN 

READ 

10 

IORD 

**r i IS 

TB4 

BOOL 

2224 

uRR 

TRA 

S. JXIT 


END 
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